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Abstract. We consider multi-messenger constraints on very heavy dark matter (VHDM) 
from recent Fermi gamma-ray and IceCube neutrino observations of isotropic background 
radiation. Fermi data on the diffuse gamma-ray background (DGB) shows a possible unex- 
plained feature at very high energies ( VHE) , which we have called the "VHE Excess" relative 
to expectations for an attenuated power law extrapolated from lower energies. We show that 
VHDM could explain this excess, and that neutrino observations will be an important tool 
for testing this scenario. More conservatively, we derive new constraints on the properties of 
VHDM for masses of lO'^-lO^'^ GeV. These generic bounds follow from cosmic energy budget 
constraints for gamma rays and neutrinos that we developed elsewhere, based on detailed 
calculations of cosmic electromagnetic cascades and also neutrino detection rates. We show 
that combining both gamma-ray and neutrino data is essential for making the constraints 
on VHDM properties both strong and robust. In the lower mass range, our constraints on 
VHDM annihilation and decay are comparable to other results; however, our constraints 
continue to much higher masses, where they become relatively stronger. 
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1 Introduction 

A new era in high-energy multi-messenger astronomy is dawning. Camma-ray and neutrino 
observations will help probe the origins of cosmic rays, the nature of high-energy sources, 
the mechanisms that power astrophysical objects, and more. In particular, these observa- 
tions may help finally reveal the particle properties of dark matter. Very heavy dark matter 
(VHDM), with mdmC^ > TeV, could produce negligible rates in nuclear scattering exper- 
iments, due to its low number density, and might be out of reach of collider production 
experiments, due to its high mass. However, there are new opportunities to search for the 
products of its annihilation or decay, due to the increasing sensitivity of gamma-ray and 
neutrino experiments. Here the high mass of VHDM may compensate its low number den- 
sity, because so much energy (2mdmC^ or rridmC^, respectively) is released per annihilation or 
decay. 

In the GeV range, the LAT instrument onboard the Fermi satellite has found various 
extragalactic gamma-ray sources, including active galactic nuclei (ACN), star-burst galaxies 
and gamma-ray bursts (CRBs) [e.g., 1-4]. The diffuse gamma-ray background (DCB) [5], 
presumed to be cosmic, was measured by LAT in the GeV range [6, 7] and was found to be 
lower than that obtained by EGRET [8]. The origin of the DGB is unsettled; in any case, 
the measured spectrum leads to strong bounds on the gamma-ray emissivity from general 
sources in the universe [9-11], as well as in specific scenarios, e.g., fast redshift-evolution 
models of ultra-high-energy (UHE) cosmic rays [12]. 

In the very-high-energy (VHE; > 0.1 TeV) range, thanks to the big successes of imaging 
atmospheric Cherenkov telescopes (lACTs) such as HESS, MAGIC and VERITAS, not only 
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Galactic but also extragalactic gamma-ray sources have been discovered [see a recent review, 
e.g., 13, and references therein]. The next generation lACT, Cherenkov Telescope Array 
(CTA), is also being planned [14]; its sensitivity is expected to be improved by a factor of 
~ 10 compared to existing lACTs. However, the discovery power of gamma-ray experiments 
at the highest energies will always be limited by the opacity of the universe caused by pair 
production processes on the extragalactic background light (EBL) and the cosmic microwave 
background (CMB). 

Neutrinos, however, are not attenuated in the cosmos, so are especially important to 
reveal distant, high-energy processes directly. The IceCube detector at the South Pole was 
completed [e.g., 15-19] and the comparable KMSNet detector is being planned in the Mediter- 
ranean Sea [20]. Interesting constraints on UHE neutrino fluxes have been placed with, e.g., 
the balloon-based ANITA experiment [21]. Although non-terrestrial high-energy neutrinos 
have not been detected so far, the sensitivity of present experiments is nearing well-motivated 
theoretical expectations [22, 23]. 

Various searches for dark matter have been made to find out indirect signals from its 
annihilation [e.g., 24-30] and decay [e.g., 31-37], including possible signatures in the DGB 
and the neutrino background. Despite the fact that dark matter has been known since the 
1930s [38] by its gravitational effects, its particle properties remain unknown. Numerous 
candidates have been suggested as a solution to this dark matter problem. Among various 
possibilities, weakly interacting massive particles (WIMPs) are especially popular [39]. They 
are motivated by the problem of electroweak symmetry breaking, and neutral WIMPs like 
the neutralino have been studied in great detail. Neutralinos are spin-1/2 Majorana fermions, 
predicted in the supersymmetric extension of the standard theory with R-parity conservation. 
They can naturally be the stable lightest supersymmetric particle, and may have annihilation 
cross sections such that their thermal relic abundance matches the observed dark matter 
density. This implies that WIMP-matter interactions may be strong enough that WIMPs will 
be produced at particle accelerators and detected directly in nuclear scattering experiments. 
Self-annihilations of WIMPs in the present-day universe may lead to detectable gamma-ray 
and neutrino signals; their non-observation has led to a variety of limits on properties of 
annihilating WIMPs [40-47]. 

Furthermore, there is no fundamental objection to considering unstable dark matter, 
as long as the lifetime is longer than the age of the universe. This possibility has been 
studied and constrained in various aspects from many years ago [e.g., 31, 32], including more 
general cases where neutral relics are subdominant [48-51]. Later, some specific particle 
physics models of decaying dark matter have recently been proposed [e.g., 34-37]. One of 
the long-lived dark matter candidates in the supersymmetric theory is the gravitino in R- 
parity-breaking vacua, where long lifetimes are allowed by the supersymmetric breaking scale 
and small R-parity violation [34, 35]. Not only the gravitino but also the sneutrino can be 
viable decaying dark matter candidates as super- WIMPs [56] that may not be seen by collider 
production and nuclear scattering experiments. It is also possible to consider hidden sector 
gauge bosons, gauginos, and fermions as decaying dark matter [36]. Composites of strongly- 
interacting particles have also been proposed [37]. Decaying dark matter scenarios have been 
constrained from the pre-Fermi era [48-53], and they have been recently discussed especially 
in the literature of the anomaly in the positron spectrum [47, 54, 55]. 

For ordinary thermal relics, the required annihilation cross section is < av >dm ~ 
(2 — 5) X 10~^^ cm^ s~^ [see 57, and references therein]. However, dark matter may be 
nonthermally produced in the early universe, and the cross section can be larger in principle. 
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Nonthermal production is commonly considered in the models of decaying dark matter. While 
many models, including supersymmetric models, consider dark matter masses of < 100 TeV, 
super heavy dark matter, e.g., wimpzillas, has also been suggested in various contexts [58-62] 
including earlier works motivated by the string theory [63]. There are various production 
mechanisms, and the wimpzilla mass, ~ 10^ — 10^^ GeV, may be related to the inflaton 
mass [60, 61]. Long lifetimes can be realized by considering breaking discrete gauge symme- 
tries and some other non-perturbative effects [62, 64]. Although the motivation to explain 
UHE cosmic rays seems to have almost disappeared [58, 62], super heavy dark matter is still 
viable dark matter candidate. 

In this work, we focus on generic candidates for present-day decaying VHDM with 
masses between 10^ and 10^° GeV. Since dark matter has not yet been discovered in the 
"expected" lower mass range, we are motivated to search at higher energies, which are now 
being explored with much greater sensitivities. There are also several independent works 
focusing on decaying VHDM [65]. Another new point of our paper is to discuss the origin 
of the possible "VHE Excess" in the DGB [10], which has been identified in the Fermi era. 
Dark matter with masses of ~ 1 — 10 GeV was discussed many years ago to explain the DGB 
in the MeV range [31]. Thanks to Fermi, the VHE DGB is now seen well above the energy 
ranges of COMPTEL and EGRET, and we find that this new feature could be explained by 
dark matter, and that IceCube searches for neutrinos are already constraining some of the 
cases in this scenario. 

Many previous studies on gamma rays focused on the mass range below 10 TeV, where 
the gamma-ray cascade caused by pair production on the cosmic photon backgrounds is 
not relevant. At higher masses, however, one has to take gamma-ray cascades into account 
properly. Our detailed calculations allow us to derive the new cascade gamma-ray bound on 
the annihilation cross section and decay lifetime of VHDM. Importantly, the Fermi DGB is 
lower than the EGRET one, so past constraints on VHDM [e.g., 48, 49, 51] can be improved 
by about a factor ~ 10. In addition, though neutrino constraints were considered in early 
days [e.g., 48-50], IceCube has indeed allow us to achieve limits that are stronger by orders 
of magnitude than those from the first generation experiments in the COMPTEL/EGRET 
era, e.g., from the Monopole Astrophysics and Cosmic Ray Observatory (MACRO) and the 
Irvine-Michigan-Brookhaven (1MB) detector. Updated calculations of both the messengers 
are of importance and enable us to demonstrate that neutrino observations can be more 
stringent than gamma-ray observations at sufficiently high masses. 

The organization of this paper is as follows. In Section 2, we revisit constraints on the 
cosmic energy budget of gamma rays and neutrinos, based on results of Fermi and IceCube, 
providing an overview of how we restrict the properties of VHDM. In Section 3, we show 
numerical results of the spectra arising from annihilating or decay VHDM. In Section 4, 
we consider the most optimistic cases, where VHDM is assumed to contribute to the VHE 
DGB, and demonstrate the importance of neutrino observations. In Section 5, we develop 
and synthesize our combined gamma-ray and neutrino constraints on annihilation and decay 
of VHDM. In Section 6, we present our conclusions. Two appendices present further details 
of the calculation of intergalactic electromagnetic cascades and how the Galactic foreground 
signals from dark matter are taken into account in our calculations of cosmic signals from 
VHDM. In this paper, according to WMAP seven-year results [66], we adopt Hq = lOO/i = 
70.2 km s"^ Mpc-\ Qdm = 0.229, Qm = 0.275 and Oa = 0.725. 
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2 Overview of cosmic energy budget considerations 



Indirect dark matter signatures from annihilation or decay can be searched for using diffuse 
background radiation. The isotropic DGB has been measured, but its origin is not under- 
stood, so a dark matter contribution can be recognized only if it is significant relative to the 
full DGB, and not merely to its uncertainty. If the DGB mainly consists of relatively rare 
point sources such as blazars, the DGB can be reduced by further progress in resolving point 
sources; the improvement is expected to be modest [43]. The neutrino background has not 
been detected so far, though atmospheric neutrinos have been measured up to several hun- 
dred TeV. At higher energies, where the atmospheric neutrino flux is negligible, the bounds 
on dark matter properties are already powerful and will improve quickly. 

It is useful to convert constraints from diffuse background fluxes into those on cosmic 
energy budgets. Gamma rays and neutrinos may be produced by astrophysical sources or 
dark matter. In the case of dark matter annihilation and decay, these stable standard model 
particles will be among the dominant endpoints of all decay chains from heavier particles. 
Both neutrinos and (lower-energy) gamma rays travel long distances, so that they are sensitive 
to processes occurring throughout the cosmos. One cannot see as far with higher-energy 
gamma rays, because they induce electromagnetic cascades, as will electrons and positrons; 
all of these are included in our calculations. 

Following Murase et al. [10], we briefly discuss the required or allowed energy budgets in 
gamma rays and neutrinos. We denote the bolometric energy budget and differential energy 
budget by Q{z) and EQe{z), respectively. For discrete sources, Q{z) is the product of the 
luminosity L{z) per source and the source density ns{z), and EQe{z) is the same before 
integration over energies; for continuum sources, these generalize. We consider two typical 
cases. For a quasi-differential limit, we consider a near-mono-energetic injection spectrum. 
For an integrated limit, we consider a E~'^ spectrum (i.e., EQe oc const.) spanning several 
decades. Constraints are reported in terms of the local {z = 0) energy budgets, with the 
assumed redshift evolution noted. 

Sufficiently high-energy gamma rays from distant sources cannot avoid pair creation 
with the EBL and CMB, and will be cascaded down to GeV-TeV energies through pair 
creation and inverse Compton scattering processes. Whatever the origin of the DGB is, if 
the cascade is sufficiently developed, it has a near-universal form [10, 67], 



where the normalization of Ge^ is / (IE^Ge., = 1- Here, E™'^ is the energy where the sup- 
pression due to pair creation occurs, /3 is typically ~ 2, and E}^^ (4/3)(£''J^"*/meC^)^ecMB — 
0.034 GeV (S^'^VO.l TeV)^((l + z)/2)^ where ecMB is the typical CMB energy. 
The DGB observed by LAT can be fitted by a power law [6], 



where = / dE^ Q e^ is the local gamma-ray energy budget and is the pre-factor coming 
from its evolution [69]. For dark matter annihilation or decay, one has ~ 1 as a typical 
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E^^^ = 0.855 X 10"^ GeV cm'^ s"^ sr'^ {E^/100 GeV)"°-^\ 
which can be crudely connected to the theoretical background gamma-ray flux as 



(2.2) 



(2.3) 
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Figure 1. Upper bounds on the local {z = 0) 
gamma-ray energy budget of the universe. The 
upper solid curves represent quasi-differential 
limits obtained with E'Qe' = Q■y^ i.e., the 
same energy input at any E' , for near-mono- 
energetic gamma-ray injection. The lower 
dot-dashed curves represent integrated limits 
for E'Qe' ~ const., i.e., an E~^ differential 
spectrum, where -E'min ~ 10^'^^ GeV and 
E'^^^ = 10^^'^^ GeV arc assumed. Thick curves 
are obtained for star formation evolution [68], 
whereas thin curves arc for no redshift evolution. 
Adapted from Ref. [10]. 



Figure 2. Upper bounds on the local {z = 0) 
neutrino energy budget of the universe. The 
upper dashed curves represent quasi-differential 
limits estimated from the IceCube-40 analy- 
sis (333.5 days) [18]. The lower dot-dashed 
curves represent integrated limits obtained for 
EuQe^ = const, (and a narrower energy range 
than that used in Figure 1). The thick curves 
arc obtained for star formation evolution [68], 
whereas the thin curves arc for no redshift evo- 
lution. The dotted curves show limits from the 
atmospheric neutrino background. WB represents 
the Waxman-Bahcall bound [69] obtained with 
EctQe^ = 0.6 X 10^" erg Mpc"^ yr"!, and MB 
represents the effective nucleus-survival bound 
of Murase & Beacom [70]. Adapted from Ref. [10]. 



value, because the most important contributions come from z ~ as long as the dark 
matter lifetime is longer than the Hubble time tfj ^. For astrophysical sources, important 
contributions often come from sources at z ~ 1 — 2, where the typical redshift of z ~ 1 
leads to E^^^l^r^i ~ 0.1 TeV, so that one has ^2 ~ a few. Since cascade gamma rays have a 
characteristic spectrum and we typically have E^Ge^ ~ 0.1, details of primary spectra are 
not very relevant because they are washed out. The DGB therefore probes the bolometric 
electromagnetic energy budget in the VHE/UHE range. From Eq. (2.3), the DGB obtained 
by LAT implies 

= L^ris < 2 X 10^5 erg Mpc-^ yr-^CS (2.4) 

when the injection energy is high enough to induce cascades. Our derived DGB constraint 
is shown in Figure 1. 

Unlike gamma rays, neutrinos reach the Earth through intergalactic space without at- 
tenuation. Because of this, neutrino detectors are directly sensitive to the differential energy 
budget and spectrum in the VHE/UHE range, EQe- Any non-terrestrial neutrino back- 
ground has not been detected yet, but the IceCube-40 integrated limit obtained for the E~'^ 

^In the no redshift evolution case that is expected for decaying dark matter, we obtain ~ 0.59. For 
annihilating dark matter, we get ~ 0.86 with the fitting function used in Yuksel et al. [46]. 
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spectrum has now reached an impressive sensitivity [18], 



El^^ < a few X 10"^ GeV cm^^ sr'^ (2.5) 

the typical differential sensitivity is weaker by a factor of a few. The theoretical extragalactic 
neutrino background is related to the energy budget as [69] 

El^u « ^{EuQE.)iz- (2.6) 
Here Ej^Qe^ is the differential neutrino energy budget, which is constrained as 

E.Qe. < 4.3 X 10^3 Mpc-3 yr-i ( fp^^ „^ ^ C'- (2.7) 

\3 X 10 GeV cm ^ s ^ sr ^/ 

The neutrino constraint is shown in Figure 2. Note that the neutrino flux limits can be below 
the gamma-ray measurements, so neutrinos are more sensitive in some cases, especially at 
high energies. At low energies, the neutrino sensitivity is worse due to the atmospheric 
neutrino background. 

It is often useful to compare the experimental neutrinos limits to theoretical bounds used 
in the literature of cosmic rays. In Figure 2, we also show the nucleon-survival (Waxman- 
Bahcall) bound derived for cosmic-ray sources [69] and the nucleus-survival (Murase-Beacom) 
bound for sources of cosmic-ray nuclei [70]. The cosmogenic neutrino flux, which gives a 
guaranteed astrophysical neutrino background, practically obeys the nucleon-survival bound 
(for primary protons) and the nucleus-survival bound (for primary nuclei) [71]. 

Our general constraints on the energy budgets can be adapted for dark matter. For 
decaying dark matter, the local energy budget, Qdm, is 

2 2 

r^dcc _ Pdra TT^dmC _ PdmC , . 

Vdm — — — — 5 l^-oj 

''''dm 'dm 'dm 

because each decay injects an energy mdmC^- For annihilating dark matter, one has 

\mdmj 2 TTldm 

because each annihilation injects an energy of 2m(jmC^- Unlike for decay, here the clustering 
of dark matter is important, and go is the flux-multiplier at z = 0, which depends on the 
dark matter profile and substructures (see Section 3.2). 



3 Calculations of neutrino and gamma-ray spectra from VHDM 

The origin of the DGB is a mystery, and a number of models have been suggested so far [see 
recent reviews 1, 5, and references therein]. Whatever the origin is, in the standard theory 
(unless axion conversion or Lorentz-invariance violation is invoked), extragalactic gamma 
rays must be attenuated at sufficiently high energies due to the pair creation with EBL pho- 
tons. Therefore, one should expect suppression at > 100 GeV energies, so that, even if the 
originally-emitted DGB extends to higher energies, it will be severely attenuated. Although 
careful data analyses and precise measurements of the VHE DGB are required, preliminary 



-6- 



Fermi data ^ suggest the possible existence of a distinct component at > 100 GeV [10], 
because the measurements are well above expectations for an attenuated power law extrap- 
olated from lower energies. We have called the discrepancy the "VHE Excess" and shown 
that it can have a significant contribution from cascade components [10]. 

The suppression due to the EBL could be compensated by an additional emission com- 
ponent if the primary spectrum of injected gamma rays is hard enough. Not including cascade 
effects, the photon index in the VHE range would need to be a < 2. However, for such a 
hard spectrum, the effects of cascades moving power from higher to lower energies cannot be 
neglected once there is enough luminosity at sufficiently high energies. This may be realized 
by some population of blazars with a very hard spectrum. Another possibility is that the 
cascades are initiated by VHE cosmic rays instead of gamma rays. It is also possible that 
the "VHE Excess" is not due to cascades, for example if the excess gamma rays are due to 
unaccounted-for foreground emission within the Milky Way, where gamma-ray attenuation 
can be ignored. 

VHDM annihilation or decay can lead to gamma-ray spectra that could explain the 
data [c.f. 73]. Extragalactic components would be cascaded, and could lead to similar final 
spectra as astrophysical sources, because cascade effects obscure the primary spectra. Fur- 
ther, there could also be a Galactic component, which would not be attenuated or cascaded, 
and this can play an important role (see Section 4.1). 

Neutrinos can reach the Earth without attenuation. Although the neutrino background 
has not been measured yet, neutrinos from dark matter annihilation or decay, as well as 
astrophysical neutrinos, may contribute to the background. The resulting neutrino spectra 
are expected to be similar to those of the emission spectra that may extend to quite high 
energies, which can in principle reveal types of allowed channels. As we stress below, taking 
into account both the gamma-ray and neutrino spectra is essential to probing VHDM models. 

3.1 Calculation of cosmic backgrounds 

Now we describe how cosmic neutrino and gamma-ray backgrounds from VHDM are calcu- 
lated. Without attenuation or cascades, injected energies are higher than observed energies 
by E' = [1 + z)E. Here we focus on the extragalactic contribution, and below we note the 
importance of the foreground component from the Milky Way. For decaying dark matter, 
the received intensity (differential flux per energy, area, time, and solid angle) is 

where dS/dE' is the primary spectrum. For annihilating dark matter, one has 

c r g{z){l + zf < av >dr,,pl^ dS 

'^""^J ^n, + \i + zfn^ dE'- 

For the decay case, the clustering of dark matter is irrelevant for the extragalactic contribu- 
tion, whereas it is very important for the annihilation case. The extragalactic clustering is 
represented by the flux-multiplier g{z), which is [74] 

g{z) = [ dM ^9{ciM,z))^^. (3.3) 

J dM Pdm "dm 



See the presentation by M. Ackermann on behalf of the Fermi coUaboration: 
http://agenda. albanova.se/conferenceDisplay. py?confld=2600. 
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Here driYie^io/ dM is the halo mass function, c{M,z) is the concentration parameter and Ac 
is the overdensity relative to the critical density pc- With the Navarro- Frenk- White (NFW) 
profile [75], for the mass function obtained by Jenkins et al. [76], go = g{z = 0) ~ 10^ is 
obtained. For the Press-Schechter mass function, 50 — (4 — 5) x 10^ [46]. However, this 
factor may be largely enhanced by substructures. For example, one of the most recent N- 
body simulations, Millenium H [77], leads to ~ 10^ — 10^ when the luminosity function is 
extrapolated down to a damping scale mass limit of 10~^/i~^A/q. This enhancement factor 
due to substructures can take a broad range of values [see also 28, 42, 78]. We simply take 
qq = 10^ and we use the fitting function given by Yuksel et al. [46] for redshift evolution. 
Other cases of substructure models can be easily compared because they scale the fluxes by 
a constant factor. 

For the neutrino background, we numerically calculate it through Eqs. (3.1) and (3.2) 
for dS/dE' obtained using PYTHIA [72]. For the extragalactic contribution to the DGB, 
we must calculate intergalactic cascades, and we perform detailed calculations by solving 
the Boltzmann equations taking into account pair creation, inverse Compton emission, syn- 
chrotron emission, and adiabatic loss (see Appendix A). Such detailed calculations are time- 
consuming but are important, especially in the high mass range, where VHE/UHE gamma 
rays are cascaded down as a result of many steps in the cascade process. Then the extra- 
galactic gamma-ray background flux is evaluated from 



where the near-universal secondary spectrum observed at the Earth is used rather than the 
primary spectrum in the cosmic rest frame. 

Finally, the resulting background fluxes can be compared to the measured DGB or the 
sensitivity curves of neutrino detectors. In the next section, we show that the expected 
theoretical gamma-ray background can indeed be compatible to the VHE DGB measured by 
Fermi. In such scenarios, the corresponding neutrino background fluxes in some cases are 
high enough to detect with IceCube within a few years. 

3.2 Choices of final states 

The final states following dark matter annihilation or decay are unknown and model-dependent 
We try to develop general limits by considering some typical final states. For the dark mat- 
ter masses considered here, all standard model particles are kinematically available, and one 
generally expects that all possibilities will occur to some degree, though there may be some 
dominant channels. As representative channels, we consider examples of gauge bosons, heavy 
quarks, and heavy leptons: (66, W~^W~, p^^~), each with 100% branching fractions. (It is 
trivial to scale our results for any other assumed branching fractions.) 

Once we show the results for these cases, we can discuss their similarities and differences, 
and say how to draw reasonably model-independent conclusions. Of course, eventual more- 
detailed studies could further these conclusions. The dependence of the spectra on mass exists 
but is not very strong, in the sense that constraints are sensitive to the typical energy flux at 
the typical energy. In the neutrino case, the differences are visible but are moderate; in the 
gamma-ray case, the differences among the assumed final states can be much weaker due to 
cascades. Thus, the choices made for the final states, while not easily justified in advance, are 
shown to be quite adequate. For ease of comparison, we apply the same particle-antiparticle 
final states used for annihilation also to decay, as expected for scalar dark matter, even 
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Figure 3. Final state spectra of gamma 
rays/clcctrons/positrons (summed) for DM — > 
DM W+W-, and DM bb with 
^/i = 100 TcV. 




log(E/s' 



Figure 4. Final state spectra of neutrinos 
(summed over flavors) for DM — fJ'~^IJ'~, 
DM W+W-, and DM ^ bb with 

^/i = 100 TeV. 



though these may not always be possible in other cases. However, the results for a single one 
of those particles, or a similar particle in that class, would be very similar. 

Some possible final states could seem quite different from the ones we assumed. Either 
77 or e~^e~ would yield no neutrino constraints, though they would have similar cascade 
bounds as each other. Seemingly, vv' would only lead to neutrino constraints but not gamma- 
ray constraints; however, electroweak bremsstrahlung leads to non-negligible electromagnetic 
branching ratios, and so gamma-ray constraints can be applied [40]. (Note that those gamma- 
ray limits based on EGRET would be improved by ~ 10 by using Fermi data.) 

To evaluate particle yields resulting from annihilation and decay, we run PYTHIA [72] . 
At > PeV energy-scales, the yields at lower energies are scaled to provide a simple approxi- 
mation to what might happen at higher energies. While this is not perfectly correct due to 
the scaling violation or new physics beyond the standard model, our results show that this 
is sufficient for our purposes. 

Examples of final state spectra are shown in Figures 3 and 4. They are consistent with 
spectra shown in previous works [e.g., 17], and satisfaction of energy conservation is checked. 
Because only a fraction of the energy goes to baryons, the amount of electromagnetic and 
neutrino energy per annihilation or decay almost matches mdmC^ in the decay case or 2mdmC^ 
in the annihilation case. For the W~^W~ channel, one sees a peak due to weak boson decay 
into a neutrino and lepton. In addition, we consider the pure leptonic channel. Such 

channels may be realized in some models such as a Sommerfeld-enhanced model including 
light force-scale carriers [25]. In Figure 3, we show the sum of contributions from gamma rays 
and electron-positron pairs, as appropriate for setting simple cascade bounds (see Section 5.2); 
in calculating secondary gamma-ray spectra below, we treat the primary gamma-ray and 
electron-positron spectra separately. 

4 VHDM signatures in the diffuse gamma-ray and neutrino backgrounds 

In this section, we show the results of gamma-ray and neutrino spectra calculated through 
the formulas provided in Section 3. In particular, we compare the theoretical background flux 
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to the measured DGB, and give some possible interpretations, first considering annihilation 
and then decay. 

Throughout the paper, the low-energy part of the DGB is assumed to be due to as- 
trophysical sources such as AGN and star-burst/star-forming galaxies. For simplicity, we 
assume that it continues to higher energies and is attenuated (its cascade effects are negligi- 
ble, because the spectrum is steep [e.g., 79]), and we adopt the following parameterization, 

^2$^ = 0.855a X 10"^ GeV cm'^ s"^ sr"^ {E^/100 GeV)^""e-^^"(-^^'"), (4.1) 

where T^^{E^,a) is the effective intergalactic optical depth that is calculated via r^^(-E^,z). 
The real situation can be complicated by having several astrophysical contributions, especially 
from blazars and star-burst galaxies [see reviews 1, 5, 79]. Complete results would depend on 
many details, luminosity functions, spectral energy distributions, and so on. To avoid such 
complexities and uncertainties, the above form is assumed to be the sum of all astrophysical 
source types. This is enough to demonstrate the potential importance of VHDM in the VHE 
DGB. Detailed modeling and precise fitting procedures are beyond the scope of this work. 

The discussions provided in this section should be regarded as the most optimistic case, 
in which the gamma-ray spectra from VHDM are a significant component of the observed 
DGB. The constraints obtained in the next section are much more general and conservative, 
in that the allowed parameter space is defined by dark matter signals being less significant. 

4.1 Possible dark matter scenarios 

We first consider annihilation of VHDM. As an example. Figure 5 shows that the annihi- 
lation of ~ 5 TeV dark matter could seemingly explain the "VHE Excess" in the DGB. In 
this case, though the attenuation is moderately important, the cascade effect is not much 
relevant because secondary gamma rays appear at < 10 GeV. The attenuation by the EBL 
is compensated by the fact that the initial spectrum is quite hard (c.f. Figure 3). Note that 
the implied cross section and the boost factor have to be quite large. Even for go = 10^, the 
corresponding cross section is < av >dm = 9.5 x lO"^^ cm^ s~^. 

As another example, Figure 6 shows a case where the intergalactic cascade is essential. 
Information on the primary spectrum is smeared out, and the results do not depend on the 
dark matter mass as long as dark matter is heavy enough. Note that, in Figures 5 and 6, 
the dark matter contributions are basically subdominant and relevant only at > 20 GeV. 
Therefore, the required cross sections are slightly below current constraints on the basis of 
the Fermi data on the DGB [42, 43]. 

In both the cases, quite large values of the boost factor and/or the annihilation cross 
section are required in order that VHDM annihilation makes a significant contribution to 
the DGB. The annihilation cross section could be larger than canonical values in some mod- 
els, e.g., the Sommerfeld-enhanced model [25], when dark matter is nonthermally produced. 
Large boost factors also seem necessary, which can be realized by significant substructure con- 
tributions. The substructure can enhance both the Galactic and extragalactic contributions, 
and the extragalactic contribution can be more important in some optimistic substructure 
models [see 42, and references therein]. For demonstrative purposes, we consider only the 
extragalactic contribution, which also gives more conservative constraints compared to the 
case where both the contributions are included. However, since the large cross section and/or 
large boost factor are required, they may conflict with other constraints [44, 80, 81]. 

Because of the above caveats, in this paper, we mainly focus on decaying VHDM. For 
decaying dark matter, the flux is insensitive to substructures, and both the Galactic and 
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Figure 5. The isotropic DGB (sum; solid 
curve) from dark matter (dashed curve) and an 
EBL-attenuated power-law extragalactic compo- 
nent (dotted curve). The assumed dark matter 
energy budget is Qdm ~ 3.0 x 10^^ erg Mpc"'^ yr""'^ 
for the W~^W~ channel. For annihilating dark 
matter, the corresponding dark matter mass is 
"T-dmC^ = 5 TeV. For the EBL-attenuated power- 
law component, a = 0.9 and a = 2.56 are assumed 
with star formation evolution. 




E [GeV] 

Figure 6. The same as Figure 5, but 
Qdm = 3.8 X 10''^ erg Mpc"^ yr'^ is used 
for the bb channel. For annihilating dark 
matter, the corresponding dark matter mass is 
'TidmC^ — 5 EcV, but resulting secondary spectra 
do not depend on the mass or on detailed primary 
spectra due to the intergalactic cascade. For the 
EBL-attenuated power-law component, a = 0.8 
and a = 2.56 are assumed. 



extragalactic contributions should be relevant. In Appendix B, we describe how the Galactic 
contributions are calculated. Introducing the J factor, i.e., the intensity of the Galactic 
contribution that is proportional to the line of sight integration of the dark matter density, 
we can define J^q as the quantity averaged over the sky region [46], which is used for the 
estimate of the Galactic contribution. 

In Figure 7, we show an example of a DGB produced by decaying VHDM with mdmC^ — 
3 TeV for the DM — t- W~^W~ channel. It suggests that the sum of the two components 
can explain the VHE DGB if the dark matter mass is in the appropriate range and the 
extragalactic astrophysical component is expressed by a simple power law ^. Though the 
Galactic component is naively expected to be comparable to the extragalactic one, it is 
actually dominant here. As indicated in Figure 7, the extragalactic component is suppressed 
by interactions with the EBL at > 100 GeV energies. On the other hand, secondary gamma 
rays by primary electrons/positrons and absorbed < TeV gamma rays are radiated in the 
GeV range, where the astrophysical power-law component is supposed to be dominant. 

We show the cases with DM — )■ /-f^/U" in Figures 8 and 9, where the Galactic component 
is dominant in the relevant energy range. In Figure 8, the hard spectral component in the 
~ 100 GeV range is dominated by secondary inverse Compton emission rather than primary 
gamma-ray emission from decaying dark matter. In Figure 9, the Galactic synchrotron 
component is dominant since UHE pairs mainly cool via synchrotron emission rather than 
inverse Compton emission. The latter situation is possible only when VHDM is super heavy. 
In both the cases, the extragalactic component also has a similar spectrum, though it has a 
tail due to the intergalactic cascade. 

Alternatively, as in Figure 10, one may consider cases where the extragalactic contri- 

^It is also consistent with present constraints from anti-proton measurements [82]. 
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Figure 7. The DGB (sum; solid curve) from 
decaying dark matter (dashed curves) and an 
EBL-attenuated power-law extragalactic compo- 
nent (dotted curve). We use JTidmC^ — 3 TeV and 
Tdm = 1-2 X 10^^ s for DM ^ W+W'. Both 
the extragalactic (thick curve) and galactic com- 
ponents (thin curve) are shown. For the power- 
law, a = 0.9 and a = 2.50 are assumed with star 
formation evolution. 



Figure 8. The same as Figure 7, but for 
mdmc^ ~ 30 TeV and Tdm = 2.0 x 10^^ s, and 
DM — > For the power-law component, 

a — 0.9 and a = 2.48 are assumed. 





E [GeV] 



10"= 
E [GeV] 



Figure 9. The same as Figure 7, but for mdmC^ = Figure 10. The same as Figure 7, but for 
10 EeV and Tdm = 1-8 x 10^^ s, and DM fi+ . nMmC^ = 10 EeV, Tdm = 4.7 x 10^^ s and DM 
For the power-law component, a = 0.9 and a = bb. For the power-law component, a = 0.8 and 



2.48 are assumed. 



a = 2.56 are assumed. 



bution is more relevant. When the Galactic component is present, this is possible only when 
dark matter is so heavy that the Galactic component exceeds the DGB measured by Fermi. 
In this case, as in Figure 6, the gamma-ray emission at ~ 100 GeV energies is dominated by 
the VHDM-induced intergalactic cascade. Since the cascade spectrum is not much sensitive 
to the primary spectrum, this intergalactic cascade scenario is allowed only in exotic cases 
where dark matter mass is higher than ~ 10 PeV. Well-motivated particle models typically 
lie below ~ 100 TeV, but super heavy dark matter is still viable and the intergalactic cascade 
spectrum is similar even for masses > 10^ GeV. 
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Figure 11. The isotropic DGBs (dashed curves) 
and neutrino backgrounds (dot-dashed curves) of 
decaying dark matter. The dark matter parame- 
ters are the same as those in Figure 7. Both the 
extragalactic (thick curves) and galactic compo- 
nents (thin curves) arc shown. The atmospheric 
neutrino background is also shown (dotted curve). 
This case is allowed. 




E [GeV] 

Figure 12. The same as Figure 11, but the dark 
matter parameters come from Figure 8. This 
case is allowed now, but is testable by the full 
IceCube. 




E [GeV] 

Figure 13. The same as Figure 11, but the dark 
matter parameters come from Figure 9. This case 
is ruled out by neutrino observations. 




E [GeV] 

Figure 14. The same as Figure 11, but the dark 
matter parameters come from Figure 10. This 
case is ruled out by neutrino observations. 



4.2 Importance of neutrino observations 

In the previous subsection, we discussed possible cases where VHDM contributes significantly 
to the VHE DGB. It would be interesting to consider how we can test those scenarios. As 
discussed in Section 2, for channels involving quarks and massive bosons, pions and other 
mesons are produced as a result of hadronization, and these lead to neutrinos. For heavy 
leptonic channels, neutrinos are produced via decay of muons and taus. Those neutrinos 
should not violate neutrino constraints as discussed in Section 3. 

In Figures 11-14, we show neutrino spectra in addition to the gamma-ray spectra (note 
the larger energy range). In the Galactic scenario demonstrated in Figures 11 and 12, the 
required dark matter mass is in the range of ~ 3 — 30 TeV. Then the corresponding neutrino 



- 13 - 



background is below the atmospheric neutrino flux, and the neutrino limits do not apply. 
In the case of DM — t- W^W" , neutrinos would not be seen by IceCube (see Section 5.2). 
The case of DM — t- may be more interesting. While it is still consistent with recent 

limits [17], the signal can be seen by the full IceCube. The non-detection can rule out this 
possibility within three years, as indicated in Section 5.2. 

As suggested in Figures 13 and 14, neutrino observations are more powerful when the 
dark matter mass is heavier than ~ 10 PeV. Such cases include the VHDM-induced cascade 
scenario demonstrated in Figures 14, where the intergalactic cascade component is dominant. 
Unless we consider channels without neutrinos (e.g., DM — )• e^e~), however, some of these 
exotic scenarios already violate the neutrino constraints. Indeed, the constraints provided in 
Section 5.1 suggests that, for such heavy dark matter masses, the neutrino constraints are 
more stringent than the cascade gamma-ray constraints. 

5 General gamma-ray and neutrino constraints on VHDM 

In this section, we show our general constraints on dark matter properties. These results are 
more conservative; the cases presented in the previous section are on the boundary of those 
being excluded. 

Cascade gamma-ray constraints are placed by requiring the calculated cascade gamma- 
ray spectrum not exceed the measured DGB data at any individual energy bin by more than 
given significance [42]. When dark matter mass lies in < TeV energy-scales, we find that 
the obtained constraints are consistent with the previous work [e.g., 43]. We extend the 
constraints to higher masses, by taking account of cosmic cascades. 

The gamma-ray constraints will be improved if more point sources contributing to the 
DGB are resolved, but the expected improvements are moderate [43]. On the other hand, 
neutrino constraints can be improved with time. Hereafter we focus on "forecasted" neutrino 
constraints to demonstrate the power of neutrino observations, assuming non-detections of 
neutrinos with full IceCube three-year observations. They are not really limits yet, and the 
neutrino background may be seen in the near future [22, 23]. 

How the neutrino constraints are placed depends on energy. We estimate the muon 
event rate by 



where fiT is the exposure and Acff(-E'jy) is the neutrino effective area. For the dark matter 
signals, the neutrino mixing is taken into account assuming 9i2 = 0.59 and 623 = 7r/4 [83]. For 
the background events, we assume the conventional atmospheric muon neutrino background 
to calculate A'bkg [22], with the assumption that any non-terrestrial neutrino background 
is not observed. The full IceCube effective area is assumed to be three times as that of 
IceCube-40 [15], and constraints are set by the criterion, Ngig/ Nsig + iVbkg < S, where 6 is 
the Gaussian significance. 

At low energies (at least in the < 300 TeV range), the atmospheric neutrino flux is 
dominant, and the constraints will improve as the square root of time. At higher energies, 
the atmospheric neutrino flux becomes negligible, and the constraints will improve linearly 
with time unless some other background is detected. 

At further high energies, the Earth becomes optically thick for neutrinos, detecting 
down-going and horizontal neutrinos is important. The above criterion would be rather con- 
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Figure 15. Neutrino (thick curves) and cascade gamma-ray (tliin curves) constraints on the anni- 
hilation cross section of VHDM. Non-detections in IceCube three-year observations and Fermi DGB 
measurements are used, respectively. Only the extragalactic contribution is taken into account, and 
constraints would be weaker for more pessimistic substructure models. The shaded region is defined 
by the more stringent bound between the neutrino and gamma-ray bounds, where for each we choose 
the most conservative limit among the selected limits for the three channels. 



servative. In order to have better constraints at sufficiently high masses, we adopt the quasi- 
differential sensitivity obtained by the analysis performed at extremely high energies [18]. 

5.1 Annihilation 

Our results for annihilating VHDM are shown in Figure 15. As in the previous section, 
we have assumed the large boost due to substructures, which could allow the extragalactic 
contribution to dominate over the Galactic contribution [see 43, and references therein]. 
Then, the gamma-ray constraints below < 10 TeV are consistent with previous works [e.g., 
43]. But, due to detailed calculations of intergalactic cascades, our results extend to much 
higher masses. Without this assumption, other constraints can give stronger results, though 
they are not focus of this work. For example, constraints from dwarf galaxies give < 10~^^ — 
10"^'^ cm'^ s~^ around 1 TeV, depending on channels [80]. The constraint is weaker as oc mdm 
as long as the energy flux limit is similar. Other limits including lACT constraints also give 
interesting limits [see a review 47, and references therein], and the recent HESS measurement 
in particular provided ^ 3 x 10~^^ cm^ s~^ at TeV for the quark-antiquark channel. Note 
that the isotropic DGB is difficult for lACTs to measure. For the neutrino constraints, 
previous limits were obtained for cases where the substructure contribution is moderate and 
the Galactic contribution is dominant. In such cases, the optimized searches for the Galactic 
halo give better constraints, and IceCube-22 results gave < 10~^^ cm^ s~^ for the fJ-^fJ-~ 
channel and < lO"^*' cm^ for the bb channel [17]. Although our limits shown in Figure 15 
are apparently better, they are more model-dependent. 

For the annihilation of point-like particles, whether dark matter is produced thermally 
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or nonthermally, the unitarity bound [see, e.g., 84] can be applied. When the annihilation 
cross section is s-wave dominated, one has [45] 



< — (w)" (3001^)"" (-) 

which may restrict the annihilating cross section of super heavy dark matter. We also show 
this unitarity bound in Figure 15. The unitarity bound can be avoided if the dark matter 
is not point-like. Note that dark matter with high cross sections in the present-day universe 
may have been produced nonthermally in earlier times. 

Beacom et al. [45] argued that neutrino limits give the most conservative bounds on 
the total annihilation cross section, because the neutrino limits are typically weaker than 
those with gamma rays. In the typically-considered mass range, that principle holds in this 
example. Interestingly, when a larger range of masses is considered, that can change. The 
present results show that the neutrino limits are stronger than the gamma-ray limits at very 
high masses. The reason is that the neutrino experiments are directly sensitive to VHE 
neutrinos, while sufficiently high-energy gamma rays are cascaded down to lower energies, 
where there are more gamma rays from other causes. Therefore, for sufficiently high dark 
matter masses, the logic about the most conservative way to set a limit on the total cross 
section would change, and the cascade gamma-ray bound could be used to put constraints 
on the total cross section. 

The shaded region is depicted as follows. In comparing these gamma and neutrino limits, 
we take the stronger of the two. Then, we conservatively take the worst limit among various 
channels. Note that the differences between the gamma-ray limits for the different channels 
are typically small, whereas this matters more in the neutrino limits. We are conservative 
about the channels that are uncertain, whereas we may use the stronger limit because both 
the messengers are reliable. 

For annihilating VHDM, we see that analytical estimates basically agree with the nu- 
merical calculations. Using ~ /emQdm (where /em is the bolometric fraction in the 
electromagnetic component), one has (for gamma rays) 



< (yv >dm < 2.9 X 10~2° cm^ s 



10 PeV 

0.113 X (\m X 10^ GeV cm^3y ( \^^\ ( f-^ 



(5.3) 



Also, using E^Qe^ ~ Qdm/T^u (where E^{dSu/dE^) = ^fsjlZy)^ one has (for neutrinos) 
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Note that, as seen in Figure 15, the constraints become weaker as dark matter is heavier. 

Generally speaking, the ratio of the Galactic contribution to the extragalactic contri- 
bution is expressed as [46] 

$G ^, BJnRscPl BJn 

^ ~ o ~ 0.1- ; 7T- — , (5.5 

$EG ctHpl^go^z (90/10^)^.' ^ ' 
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where B is the boost factor for the Galactic contribution. Note that, for gamma rays, the de- 
nominator can be affected by propagation effects. While the boost toward the Galactic center 
is order of unity, the overall Galactic boost factor can take a wide range of values [see, e.g., 
28], so that the above ratio highly depends on substructure models. In the previous section, 
we discuss possibilities that VHDM contributes to the DGB, where significant boosts due to 
substructures are assumed for annihilating dark matter. To avoid complexities due to Galac- 
tic substructures, for obtaining the constraints, we do not include the Galactic contribution 
in the present study, which means Figure 15 gives conservative constraints for a given qq. 
In the absence of substructures, the extragalactic contribution is typically subdominant [46], 
and the constraints from regions around the Galactic center should be more stringent [24]. 

5.2 Decay 

In Figure 16, we show our numerical results for present-day decaying VHDM. The gamma- 
ray constraints give Tdm ^ 10^^ — 10^^ s. For the DM — )• channel, our gamma-ray 
constraint is consistent with the previous result that is obtained only at < 10 TeV [e.g., 
43, 54]. Our result is applicable even at higher masses, thanks to detailed calculations 
of intergalactic cascades. One sees that the cascade gamma-ray bound is largely mass- 
independent at sufficiently high masses. This is because the energy injection rate of decaying 
dark matter does not depend on rridm (because the number density scales as oc l/mdm), and 
the cascade spectrum is near-universal so the DGB constrains the bolometric electromagnetic 
energy budget without depending on details of final state spectra. Also, the cascade bound 
is reasonably channel-independent. While this is obvious when the energy fraction into 
electromagnetic components of decay products is similar, it may not be true e.g., for the z^P 
channel even if all decay products are standard model particles. However, even in such cases, 
the cascade gamma-ray bound can be effectively applied when it can be regarded as a more 
conservative constraint compared to the neutrino constraint. 

The gamma-ray constraints are more stringent than neutrino ones at lower masses, 
whereas the neutrino constraints become stronger for masses of > 10 — 100 TeV (as above 
for Ref. [45] and annihilation, this reverses the logic for Ref. [53] and decay). In other words, 
neutrino observations give conservative limits at < 15 TeV for DM — )• /x"*"//", < 25 TeV 
for DM — >• W~^W~ , < 100 TeV for DM bb, respectively. Compared to the gamma-ray 
constraints, the neutrino constraints are sensitive to the choice of decay channels. This is 
because the differential spectrum is relevant for neutrino observations, and harder neutrino 
emission is more favored to avoid the atmospheric neutrino background and to be detected 
with a larger effective area. Our results are basically consistent with previous limits that are 
obtained at < 100 TeV where the atmospheric neutrino background is measured [53, 55], but 
we also cover the higher mass range. And our constraints are expected to be slightly better 
than those of Ref. [53] because the ratio of the signal to the background increases with time. 
Also, the forecasted limits are expected to be much better than the previous IceCube-22 
results unless the diffuse neutrino background is discovered [17]. Note that, of course, the 
lifetime for the cases shown in Figures 5-14 is consistent with the limits shown in Figure 16. 

For decaying dark matter at higher masses, the cascade gamma-ray bound becomes 
more conservative than the neutrino limits. As discussed in the previous section, the cascade 
leads to a near-universal spectrum, so the DGB constraints are sensitive to the bolometric 
energy budget rather than the differential energy budget, and details of the primary spectra 
do not matter for the constraints. Therefore, as in annihilating dark matter, the cascade 
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Figure 16. Neutrino (thick curves) and cascade gamma-ray (thin curves) constraints on the life- 
time of decaying VHDM (90 % CL). Non-detections in IceCube three- year observations and Fermi 
DGB measurements are used, respectively. Both of the Galactic and extragalactic contributions are 
included. The shaded region is defined by the more stringent bound between the neutrino and gamma- 
ray bounds, where for each we choose the most conservative limit among the selected limits for the 
three channels. 



gamma-ray bound could be used to put constraints on the total decay rate at sufficiently 
high masses. 

Our numerical results roughly agree with analytical expectations for the cosmic limits. 
Using Q-y ~ /emQdm (where /em is the bolometric fraction in the electromagnetic component), 
we have (for gamma rays) 



1 < 27 -1 ( 0-113 \ / 1.05X105 GeV cm'^ X f f^^^\ 



(5.6) 



where the constraints are expected to be largely mass-independent for present-day decaying 
VHDM [e.g., 49]. Also, using 

EuQe, ~ Qdm/'Tiu (where El{dS^/dE^) = yfs/U^,), one has 

(for neutrinos) 
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The cascade gamma-ray bound is simply understood through Eq. (5.6). For demonstra- 
tion, based on the result in Figure 1, we also apply the energy budget constraint in the no 
redshift evolution case, ~ 2.6 x 10^^ erg Mpc"^ yf-^) to decaying VHDM. Results are 
shown in Figure 17. One sees that the simple cascade gamma-ray constraints are comparable 
to the detailed constraints shown in Figure 16 especially in the ~ 100 PeV-1 EeV range. 
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Figure 17. The simple cascade gamma-ray bound on tlie lifetime of decaying VHDM. The limits 
are set by the energy budget constraints shown in Figure 1. The lack of information in this figure 
is exactly the point: the results are insensitive to dark matter mass and decay channels we consider 
here. Note that only the extragalactic contribution is taken into account. They are consistent with 
the detailed constraints shown in Figure 16. 



In Figure 16, both of the extragalactic and Galactic contributions are taken into account. 
The ratio of the Galactic component to the extragalactic component is 

— — — ~ 0.5—. (5.8) 

For conventional dark matter density profiles such as the NFW profile and Einastro pro- 
file [85], each contribution is typically comparable, though the situation would be affected if 
the dark matter is more centerally-clustered by e.g., dark matter contraction due to baryon 
dissipation [86]. In fact, in Figure 16, the Galactic component comparably contributes to the 
neutrino constraints. On the other hand, for the gamma-ray constraints at sufficiently high 
masses, only the extragalactic component matters since most of the gamma rays do not lie 
in the LAT range. However, when the mass is in the TeV range, the Galactic component 
in the VHE DGB is more important since the extragalactic component is affected by the 
EBL attenuation. The Galactic contribution is also enhanced when dark matter is super 
heavy (see Figure 16) when UHE electrons /positrons can emit GeV gamma rays via the 
synchrotron radiation. Getting the Galactic part would be relevant for evaluating the dark 
matter parameters. Indeed, the case with the Galactic contribution suggests the lower mass, 
compared to the case without it. 

In Figures 11-14, we showed neutrino spectra as well as gamma-ray spectra that are 
compatible to the VHE DGB measured by Fermi. As remarked in Section 4.3, the forecasted 
three-year neutrino constraint will be consistent with the case of DM — )• W~^W~ shown in 
Figure 11. However, the case of DM — t- fJ-^fJ-" shown in Figure 12 may contradict the neutrino 
constraint, that is, this case will be tested in a few years. The cases shown in Figures 13 and 
14 are exotic, and they are already excluded by the neutrino observations. 



- 19 - 




E [GeV] log(m [GeV]) 



Figure 18. The same as Figure 11, but the 
final state spectra are obtained from DGLAP 
equations. All the used parameters are the 
same as those in Figure 14, which indicate 
that the resulting constraints will be similar to 
those for the bb channel. For the dark matter 
component, TOdmC^ = 10 EeV, Tdm = 4.7 x 10^^ s 
and DM bb are adopted. For the power-law 
component, a ~ 0.8 and a = 2.56 are assumed. 
This case is ruled out by neutrino observations. 



Figure 19. Influences of the final state spectra 
obtained from DGLAP equations on neutrino 
(upper thick curves) and cascade gamma-ray 
(lower thin curves) constraints. Decaying VHDM 
is considered at rridmC^ > 10 PeV for DM — > qq. 
Note that the mass range is different from 
Figure 16 since we here focus on even higher 
masses. 



So far, we have considered final state spectra that are obtained by PYTHIA with extrap- 
olation to higher energies. For channels with accompanying fragmentation and hadronization, 
the final state spectra at higher energies become softer at small values of E/^/s (i.e., in the 
low-x regime) because of the scaling violation. In addition, including supersymmetry (SUSY) 
leads to much softer spectra [87-89]. In order to demonstrate that our essential results at 
y/s > 10 PeV are not much changed by such details of the spectra, we show an example 
for final state spectra obtained by Ref. [90] (that agrees with Ref. [89] within uncertainties), 
where Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations with supersymmetry 
(SUSY) are solved [c.f. 87, 91, 92]. The results are shown in Figure 18, where the spec- 
tral shape was calculated for an extreme energy scale m^t^^c^ = 10^^ GeV are used (so the 
realistic spectra would be slightly harder). The primary spectra in Figure 18 are indeed 
softer than those in Figure 14, but one sees that all the used parameters are the same as 
those in Figure 14. We also compare multi-messenger constraints on the lifetime of decaying 
VHDM between the above two cases. As shown in Figure 19, the constraints are similar 
within about a factor of two. The reason is that the constraints are not very sensitive to the 
detailed spectra at small values of E/^/s. For the neutrino constraints, in the UHE range, 
the limits are mostly sensitive to the energy flux around the peak, which does not vary much 
among various models. The firmness is enhanced for the gamma-ray constraints. The cascade 
gamma-ray bound is connected to the bolometric electromagnetic energy budget, since the 
initial spectral information is erased and the near-universal spectrum is achieved for suffi- 
ciently high masses. Therefore, only the energy fraction into electromagnetic components of 
decay products does matter and details of the final state spectrum are irrelevant. The energy 
fraction into electromagnetic components is quite similar among the models [92], so our re- 
sults shown in Figures 15 and 16 give reasonably accurate results within larger astrophysical 
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uncertainties. 

Note that, though we focus on present-day decaying dark matter with Tdm > tn, there 
are cosmological bounds, which are more robust but much weaker. By evahiating how the 
expansion rate is changed, one can obtain bounds on the hfetime, > 10^^ s, using CMB 
data [93]. Though the only assumption is that dark matter particles decay into relativistic 
particles, it is difficult to improve by further observations, since the decay affects CMB only 
at large scales and errors are limited by cosmic variance. Nevertheless, the analyses combined 
with type la supernovae data improve the limit by about an order of magnitude [93]. Also, 
because of our interests, the mass abundance of relics is fixed to be all of the dark matter. 
Allowing for possibilities that Tdm ^ tn and/or decaying particles are subdominant, one can 
obtain more general constraints on relic neutrals from the Big-Bang nucleosynthesis as well 
as the diffuse backgrounds [e.g., 48, 50, 51, 94]. 

6 Summary and discussion 

In this paper, we investigate diffuse neutrino and gamma-ray backgrounds coming from 
VHDM with ?TidmC^ ^ TeV, on the basis of observations by Fermi and IceCube. Our results 
are summarized as follows. 

(1) We pointed out the possibility that VHDM could contribute to the VHE DGB. This 
result is based on the recent indication of the "VHE Excess" that is identified in the Fermi 
era from the fact that the simple extrapolation of the astrophysical power-law component 
falls well below the measured DGB due to the severe EBL attenuation. We mainly considered 
decaying VHDM as the possible scenario, since quite large values of the cross section and the 
boost factor are required for annihilating VHDM. 

Although the precise determination of the VHE DGB is obviously required and other 
interpretations are also possible, significant impacts on understanding dark matter and the 
DGB are expected if the VHE Excess is confirmed. Such the optimistic scenario is also 
interesting because they are testable for some parameters. In the Galactic scenario, dark 
matter masses of mdmC^ ~ 3 — 30 TeV are necessary to explain the VHE DGB. Although the 
neutrino background is below the atmospheric background, the case of leptonic channels is 
testable by IceCube within a few years. On the other hand, in the VHDM-induced cascade 
scenario, extremely heavy masses (> 10 PeV) are required. Since large neutrino fluxes are 
predicted, such exotic cases are already excluded by neutrino observations. Those results 
demonstrate the importance of neutrino detectors such as IceCube and KM3Net for testing 
the properties of dark matter. Of course, deeper observations by Fermi and future detectors 
such as GAMMA-400 [95] should also provide us with crucial information. 

Although we considered the isotropic backgrounds for conservative studies, future more 
detailed analyses using the angular information can give tighter constraints. The dipole 
anisotropy can be expected in the Galactic scenario, so searches for anisotropy in the VHE 
range would also be useful. 

In addition, we also expect that not only neutrinos but also TeV gamma-ray observations 
via Cherenkov telescopes such as HAWC and CTA will be helpful. If high-energy gamma 
rays are produced above TeV energies due to annihilation or decay of VHDM, there may be 
a bump around 1 TeV or higher energies in the gamma-ray spectrum of galaxies. Predicted 
fluxes depend on dark matter parameters, but future searches for TeV gamma rays from 
nearby galaxies such as LMC and M31 will be relevant. 
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(2) We revisited the cosmic energy budget argument in light of Fermi and IceCube 
observations, and apphed it to VHDM. In the high mass range, gamma-ray cascades must 
be taken into account for the extragalactic contribution. We performed detailed calculations 
of both neutrinos and gamma rays, which enables us to demonstrate the power of neutrino 
observations compared to gamma-ray observations. 

We derived a new cascade gamma-ray bound both analytically and numerically. In the 
sufficiently high-mass range, this cascade gamma-ray bound gives conservative constraints, 
so that it can be used reasonably independently of channels and it is also useful to restrict 
the total decay rate and the total annihilation cross section of VHDM in the high mass 
range. Importantly, since the Fermi DGB is lower than the EGRET one, the constraints 
are improved by about a factor ~ 10. In particular, for decaying VHDM, we obtained 
Tdm ^ 10^^ — 10^'' s, which is more stringent than the neutrino limits at < 10 — 100 TeV. On 
the other hand, at higher mass scales, the neutrino limits are stronger, while the DGB gives 
more conservative constraints. The cascade gamma-ray bound is largely mass-independent at 
sufficiently high masses. For annihilating VHDM, the gamma-ray limits are typically stronger 
in the relevant mass range where the unitarity bound is preserved. But neutrino observations 
give more conservative but still interesting constraints compared to the unitarity bound. 

We focused on gamma-ray and neutrino constraints in light of Fermi and IceCube/KMSNet, 
but cosmic-ray experiments including observations of UHE gamma rays and the arrival dis- 
tribution in UHE cosmic rays are also useful in the wimpzillas regime. The importance of 
such UHE particle observations was studied in early days, and top down models for UHE 
cosmic rays have been constrained. 



Acknowledgment s 



K. M. is supported by JSPS and CCAPP. The research of J. F. B. is supported by NSF 
Grant PHY-1101216. We thank Shunsaku Horiuchi, Boaz Katz, Alexander Kusenko, Ranjan 
Laha, and especially Carsten Rott for helpful discussions. 



A Intergalactic cascades 



Taking into account the pair creation, inverse Compton, synchrotron radiation and adiabatic 
loss, we numerically calculate the cascade emission by solving the Boltzmann equations that 
are often referred as kinetic equations [87, 96, 97], 



dx 
ON, 
dx 



-N'yRjj H h 



dNr 

dx 



dx 

NeRlC + 



dx 
dx 



(A.l) 
(A.2) 
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Figure 20. Interaction length (solid curve) of 
high-energy photons for the pair creation process, 
the energy loss length of electron-positron pairs 
for the inverse Compton process (dashed curve), 
and the synchrotron cooling length for 10^'^ nG 
that can be expected in structured regions of the 
universe (dotted curve). Thick/thin curves repre- 
sent cases without/with the CRB. 



Figure 21. The average of the line of sight 
integration, inside a cone with half angle ijj 
around the Galactic center, as a function of the 
pointing angle with respect to the Galactic center 
direction for the NFW profile. The solid curve 
is for decaying dark matter, whereas the dashed 
curve is for annihilating dark matter. 



where 



Ric 

dx 

dNr 

dx 
dx 



dn f dn ^ 



dn f dn ^daic, ^ 



Here c = (1 — //)c, Pgyn is the synchrotron energy loss rate, Pad is the adiabatic energy loss 
rate, A'^^ and N^. are photon and electron/positron number densities per energy decade, and 
Q™-' and Q"^'' are photon and electron/positron injection rate. 

For the intergalactic cascade, x corresponds to travel distance (divided by c), whereas 
time is used as x for the source cascade. This paper focuses on the former case, so we 
consider the CMB and EBL get photon spectrum, dn/de. As an EBL model, we use 

the low-IR model of Ref. [98], which gives relatively conservative estimates. For variations 
in EBL models, we refer to Ref. [99]. The highly uncertain cosmic radio background is not 
taken into account in this work, but it is relevant for direct detections of UHE gamma rays 
themselves [97, 100]. Relevant length-scales are depicted in Figure 20. 
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B Galactic contribution 



In this section, we describe how to calculate the Galactic contribution of dark matter. As 
for spectra of primary particles, the Galactic contribution is given by 

rdec. ,.^ _ 1 /■ ^, Pdm{r) dS _ 1 dS 



4vr J mdmTdm dE 47rmdmTdm dE Jq 
- '^^ J'^'^'i^l:) (B.l) 



4vrmdmTdm dE 
for decaying dark matter, and 



1 f M < (^'^ > dm f Pdm{r)Y dS <av>d^dS 



U^) = — dl ^^^^^^ ^ = : — / dl pi^(r) 

< CrV >dm^scPsc dS 



dm 



STrml^ dE 



(B.2) 



for annihilating dark matter. 

We have introduced the J factor. For decaying dark matter, one has 

1 /"^max 

jdcc(^) = / dl pUr). (B.3) 

For annihilating dark matter, one has 

J-n(^) = / dl pUr) (B.4) 

^scPsc JO 



Here r = R'^^ — 2lRsc cos il^ + P, /max = -Rsc cos ^ + \/ -R^w ~ -^sc sin^ V') a-iid ^mw = 20 kpc 
is the size of our Galaxy. We take Rsc = 8.5 kpc and pscC^ = 0.3 GeV cm^'^ with the NFW 
profile. Note that J^{iIj) never vanishes, and the minimal isotropic component is written as 
^"(180°). By evaluating Eqs. (B.3) and (B.4), we obtain ~ 1 for the decay and ~ 0.4 for 
the annihilation [46], respectively. It is often useful to see the average of ^ in a cone with 
half- angle ip around the Galactic center, which is defined as 



2tt 



Jn = ^ d(cosV') W) (B.5) 

where = 27r(l — cos^) is the solid angle of a field of view, and J^n is shown in Figure 21. 

For the Galactic component, the mean free path of gamma rays is > 10 kpc, so that 
the primary gamma-ray-induced cascade would not be much critical. On the other hand, 
electrons and positirons are largely confined and should lose their energies as radiation, so that 
resulting secondary emissions have to be considered. Instead of performing time-dependent 
calculations, in this work, we obtain the stead-state distribution. In the continuous energy- 
loss approximation, the kinetic equation becomes 

8N 8 ■ ■ 

where Pic is the inverse Compton energy loss rate. Hence, the steady state distribution is 
easily obtained by 

= P Ip I Q'^' (^-7) 

-TIC + Jsyn J E 
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For a given electron/positron distribution, the gamma-ray spectrum coming from synchrotron 
and inverse-Compton emission is evaluated simultaneously from Eq. (A. 3). We adopt 1 ^uG 
as the Galactic magnetic field. 
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